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Summary 

The mean curvature of the influence surface of the space-time point ( x , t ) 
appears in linear supersonic propeller noise theory and in the Kirchhoff for- 
mula for a supersonic surface. Both these problems are governed by the linear 
wave equation with sources on a moving surface. The influence surface is 
also called the £ - surface in the aeroacoustic literature. This surface is the 
locus, in a frame fixed to the quiescent medium, of all the points of a radiat- 
ing surface /( x, t) = 0 whose acoustic signals arrive simultaneously to an 
observer at position x and at the time t. Mathematically, the £ - surface is 
produced by the intersection of the characteristic conoid of the space-time 
point ( x , t) and the moving surface. In this paper, we derive the expression 
for the local mean curvature of the £- surface of the space-time point ( x , t ) 
for a moving rigid or deformable surface /( x, t) = 0 . This expression is a 
complicated function of the geometric and kinematic parameters of the sur- 
face /( x, t) = 0. Using the results of this paper, the solution of the govern- 
ing wave equation of high speed propeller noise radiation as well as the 
Kirchhoff formula for a supersonic surface can be written as very compact 
analytic expressions. 
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1.0 Introduction 


The prediction of the noise high speed helicopter rotors and propellers is an 
important aeroacoustic problem. Until recently there have been two 
approaches in wide practice to attack this problem. The first is based on the 
Ffowcs Williams-Hawkings (FW-H) equation [1] and the second is based on 
the Kirchhoff method. The Kirchhoff method was suggested first by Hawk- 
ings for aeroacoustic applications [2]. There are now a large number of publi- 
cations on the use of FW-H equation. In general, these applications use this 
equation with surface sources specified on an impenetrable surface such as a 
blade surface. See the many references in various chapters of Hubbard [3] 
and Goldstein [4]. A good review of the applications of the Kirchhoff method 
is by Lyrintzis [5]. Ffowcs Williams also proposed to use FW-H equation 
with the data specified on a penetrable surface [6, Chap. 11, Sec. 10]. This 
was carried out for high speed helicopter rotor noise prediction by Di 
Francescantonio [7] and, independently by Brentner and Farassat [8]. This 
new method essentially brings the two approaches together and appears to be 
superior to the earlier approaches [8]. The point of relevance to the present 
paper is that in all these methods, one has to solve the linear wave equation 
with sources on a moving surface. We will address the difficult problem of 
the solution of the linear wave equation with supersonic surface sources here. 
As in all of our previous works, we will solve this problem in the time 
domain. 

The mean curvature H F of the influence surface of the space-time point 

(x, t ) appears in linear supersonic propeller noise theory and the Kirchhoff 
formula for a supersonically moving surface [9-1 1]. We refer to the influence 
surface as the E - surface in this paper to be consistent with our earlier publi- 
cations on the subject of linear wave propagation. Assuming that the propel- 
ler blade or the Kirchhoff surface is given by the moving surface 
f(x, t ) = 0, the E- surface is constructed by the intersection of the charac- 
teristic conoid of the space-time point ( x , t ) and the moving surface / = 0 . 
This construction is explained in the present paper in detail. The E - surface 
is, thus, the locus, in the frame fixed to the quiescent medium, of all the 
points on the surface / = 0 whose acoustic signals arrive simultaneously at 
x at the observer time t. Therefore, the E - surface describes the domain of 
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dependence of the space-time point ( x , t ) in the frame fixed to the quiescent 
medium. 

In our work on linear supersonic propeller noise prediction in the time 
domain, many curvature terms of the blade surface appear in the theoretical 
formulation [9-11]. Our experience with numerical evaluation of propeller 
noise indicates that the curvature terms can contribute significantly to the pre- 
dicted noise. For this reason, as well as to obtain the Kirchhoff formula for a 
supersonically moving surface, we give the detailed derivation of the analytic 
expression for the mean curvature of the £- surface. The derivation requires 
the knowledge of classical differential geometry and tensor analysis [12-14]. 
In practice, the geometry and the kinematics of the moving surface / = 0 
are known. The mean curvature H F of the £ - surface must be described in 

terms of the known geometric and kinematic parameters of the surface 
/ = 0 . This is achieved by first defining a local time dependent coordinate 
system on the £ - surface in terms of the local coordinate system on the sur- 
face / = 0 . The coefficients of the first and the second fundamental forms of 
the £ - surface are then written in terms of the known local parameters of the 
surface / = 0 . From these coefficients we construct the expression for the 
mean curvature H F . 

In section 2, we explain how the mean curvature of the £ - surface appears in 
supersonic propeller noise theory and the Kirchhoff formula for moving sur- 
faces. In section 3 we describe the construction of the £- surface and how to 
visualize this surface for a given observer space-time point (jc, t ) when the 
motion and the geometry of / = 0 are specified. In Section 4, which is the 
main part of this paper, we give the detailed derivation of the mean curvature 
H F assuming that the surface / = 0 is rigid. In a subsection, we present a 

brief derivation of H F for a deformable surface. Discussion of the main 

results and concluding remarks are in Section 5. Here we discuss the validity 
of the expression for the mean curvature and the possibility of singularity for- 
mation in the acoustic field of a supersonically moving surface. 

The main mathematical reference for this paper is [15] where all the rules of 
manipulation of multidimensional generalized functions used here are pre- 
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sented. More details on differential geometry and Green’s function solution 
of wave equation with sources on a moving surface are presented in a NASA 
Technical Memorandum [16]. We use both vector and tensor notations in the 
present paper where appropriate. The summation convention is also used 
here. 


2.0 How the Mean Curvature of the Z -Surface Appears in 

Applications 

In this section, we briefly explain how the mean curvature H F appears in two 

problems of acoustics. The first is the supersonic propeller noise prediction 
problem and the second is the derivation of the Kirchhoff formula for a super- 
sonically moving surface. Both these problems are governed by the linear 
wave equation with sources on a moving surface. 

2.1 The Supersonic Propeller Noise Problem 

Consider a supersonic propeller whose blades are described by 
f(x,t ) = 0 ,f>0 outside the body. We assume that f(x,t) = 0 is so 
defined that V/ = n , where n is the local unit outward normal to this sur- 
face. This can always be accomplished. Using the acoustic analogy of Light- 
hill, the governing equation is the Ffowcs Williams-Hawkings (FW-H) 
equation. For the purpose of discussion here, we are only interested in the 
surface source terms of this equation. The governing linear wave equation is, 
thus, the following: 

°V= !^Po v „ 6 (/)]-g^[pH,-S(/)] 

2 

where p' = p Q c , p () and c are the density and speed of sound in the undis- 
turbed medium, respectively, v is the local normal velocity on the blade and 

p is the (gage) pressure on the blade. The Dirac delta function is denoted 
b(f ) . We have assumed inviscid fluid medium here and have neglected the 
viscous forces on the fluid at the blade surface. 
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In practice, one usually divides the propeller blade surface into panels and 
then sums up the noise produced by each panel [11]. For panels where the 
Doppler factor 1 - M r is not small, a noise prediction formula based on the 

Green’s function solution of Eq.(l) with the Doppler factor (e.g., Formulation 
1A of Langley [17,18]) can be used. Here M r is the local Mach number of 

the surface in the radiation direction based on the speed of sound in the undis- 
turbed medium. The Doppler factor appears in the denominator of this solu- 
tion of Eq.(l) and, therefore, it is also referred to as the Doppler singularity. 
For supersonic propellers, the Doppler factor can be small or zero for many 
panels on the blade for a given observer position and time. In that case, 
another solution of the FW-H equation without the Doppler singularity must 
be used for noise prediction [10,11]. We present a sketch of the derivation of 
this solution below. 

We assume that the blade surface / = 0 is smooth. This is a simplifying 
assumption which is not satisfied in practice, e.g., at the trailing edge of the 
blade, but does not concern us here. See references [10,11,16] for the general 
case. Using the rules of generalized differentiation [15,16,19,20], the deriva- 
tives on the right side Eq.(l) are taken explicitly giving the following result: 

A[p 0 V„S(/)]= P 0 V„5(/)-P 0 V b 2 §'(/), (EQ2, 

A[p„.5(/)] = - 2H f pd(f) + p8(f) (EQ3) 

where Hj is the local mean curvature of the blade surface / = 0 and a tilde 

under a variable, such as p , denotes the restriction of the variable to the sur- 
face f- 0. See references [15,16] for the concept of restriction. As mentioned 
in references [9,15,16], the purpose of restriction of a variable to / = 0 is to 
reduce algebraic manipulations as well as obtaining analytic expressions in 
the simplest possible form. Using Eqs.(2) and (3) in Eq.(l), we get: 

□V= (P 0 E„ + 2H f p)8(f ) - (p + p 0 7 )8 '(/) (EQ4) 

We will give the solution of the wave equation with these types of inhomoge- 
neous source terms after the following discussion on the Kirchhoff formula 
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for a supersonically moving surface. We mention here that it is the second 
term on the right of Eq.(4) involving 8'(/) which will give an integrand 
depending on the mean curvature of the E- surface when the Green’s func- 
tion solution of Eq.(4) is sought. 


2.2 The Kirchhoff Formula For a Supersonically Moving Surface 

The need for the development of the Kirchhoff formula for a supersonically 
moving surface appears in the prediction of the high speed rotor noise of a 
helicopter. Even though the conventional rotor blades travel at subsonic 
speed, the blade tip shock can move at supersonic speed. One possible 
method of prediction of high speed rotor noise is using the Kirchhoff formula 
for a moving surface. We will call the surface on which the data are specified 
the Kirchhoff surface. In practice one specifies the data on the Kirchhoff sur- 
face from high resolution unsteady aerodynamic calculations. Kirchhoff for- 
mula based methods are showing great promise in aeroacoustics [5]. 


Let the supersonically moving Kirchhoff surface be described by the equation 
/( x, t) = 0 ,f>0 outside the surface and |V/| = n . We assume again that 
/ = 0 is a smooth surface. In practice, this assumption is not satisfied 
because it is often more convenient to use a Kirchhoff surface with edges. 
References [9] and [16] give the general result which applies to piecewise 
smooth surface. Let p' stand for the acoustic pressure satisfying the linear 
wave equation in the region outside the surface / = 0 . The governing equa- 
tion for deriving the Kirchhoff formula is obtained by first extending p' to the 
entire space as follows: 


J p /> 0 (outside the surface) 
1 0 / < 0 (inside the surface) 


(EQ5) 


Then, p' has a discontinuity across the surface / = 0. It is seen that in the 
entire three dimensional space, we have 

□ 2 /3' = 0 (EQ 6) 
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where all the derivatives in the wave operator are ordinary derivatives. The 
generalized D’Alembertian of p' gives the following equation 

[9.15.16.19.20] : 

□ 2 P' = - (p n + - c M n p\(f)-l^lM n p5(f)]-V.[pn5(f)] (EQ7) 

where the bar over the wave operator stands for generalized differentiation 

[15.19.20] . Note that in this paper we use a bar over a differentiation operator 
to signify generalized differentiation only if there is the possibility of confu- 
sion between the ordinary and the generalized derivatives of a function. For 
example, in Eqs.(l) and (4) we have not used a bar over the wave operator 
even though the derivatives in the wave operator are generalized derivatives. 
It is clear from the right sides of these equations that the wave operator can 
only be generalized D’Alembertian. 

Taking the derivatives on the right side of the above equation explicitly, and 
using results similar to those of Eqs.(2) and (3), give the following governing 
equation for deriving the Kirchhoff formula for a supersonically moving sur- 
face: 


n 2 p' = - 


Pn + \ M „P, + -M„P\ ~ 2H /Pl 8 </> 


-a -¥„ 2 ) p 5 '(/) 


(EQ 8) 


Here, as before, we have used a tilde under a variable to indicate the restric- 
tion of the variable to the surface / = 0 . Note that we can not drop the tildes 
in (M p) , where the subscript t stands for partial differentiation with 

respect to time, since ( M n P ) ^ ( M n p (see reference [15]). As in the case 

of Eq.(4), it is the second term on the right of Eq.(8) involving 8'(/) which 
gives an integrand depending on the mean curvature of the £ - surface when 
the Green’s function technique is applied to Eq.(8) to find the Kirchhoff for- 
mula for a supersonically moving surface. 

Below we will show how one interprets the solution of the wave equation 
with sources involving 8(/) and 5 ’(/) . 
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2.3 The Solution of Wave Equation With Sources Involving 8 (/) and 

«'(/) 

We now give the solution to the following two wave equations [15,20]: 

□ 2( L = q{x, t)b(f) (EQ 9) 

□ 2 <h = q(x, t)b'(f) (EQ 10) 


Let F(y; x,t) = f\y,t-~ = [f(y, x)l , , where (x, t ) and (y, x) are 


ret 

the observer and the source space-time variables, respectively. The subscript 
ret stands for the retarded time and r = \x - y| . The surface F = 0 is the 
E - surface of which more will be said in the next section. The solution of 
Eq.(9) is: 


f - 

[ q(y, T)! 

J r 

F = 0 

A 


dL 


(EQ 11) 


-I ret 


where the parameter A is defined below: 


A = (1 +M n - 2M n cos 0 ) 


1/2 


(EQ 12) 


Here, M is the local normal Mach number of the surface f - 0 and 0 is 
the angle between n and the radiation direction r = x - y. 

The solution of Eq.(10) is given by: 


4n<&(x, t ) = 


■ J 


F = 0 


dN 


+- 


if q(y,*) 

) ■ 

>s 1 

> 

Jret. 

/ q(y, t) ) 


l A J 

ret 


d£ 


(EQ 13) 


where d/dN = N.V and N is the unit outward normal to the E - surface. 
We will give this unit normal explicitly in section 4. The mean curvature of 
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the £- surface, H F , appears in the above equation. Therefore, it appears in 

solutions of Eqs.(4) and (8). Note that the mean curvature H j of the surface 

/ = 0 also appears in Eqs.(4) and (8). This curvature depends on geometry 
of the surface /= 0 which is known. We will not write the solutions to equa- 
tions (4) and (8) since we have given the complete solutions in a more general 
setting in our earlier papers [9,10,15,16]. However, note that the solutions of 
Eqs.(4) and (8) as given here are very simple and compact analytic expres- 
sions in contrast to what we have published earlier, e.g., in [11]. We will say 
more on what £ - surface is and how to visualize it. 

3.0 How To Construct and Visualize the £ - Surface 

As mentioned above the mathematical representation of the £ - surface is 
given by the relation: 

F(y;x,t) = /(y, = U(y^)] re t = 0 (EQi4) 

We give some facts and information about this surface and the way to con- 
struct and visualize it in space. The characteristic cone of the wave operator: 

, — ,2 1 d 2 3 2 

D = ^TT-Tir <EQ15) 

c di ay . 

at the point (x,t) is given by the relation: 

lx — vl r 

g = T - t + — =T-f + - = 0 (EQ 16) 

C C 

Note that the variables of the wave operator are (y, t) which are also the 
variables of the equation of the characteristic cone. These variables are called 
the source (or the field) space-time variables. The vertex of the cone is at the 
point (x, t) which is kept fixed in the analysis below. The space-time vari- 
ables (x, t ) are known as the observer variables. Equation (16) describes a 
cone in space-time because if the point (y, t) is on the surface g = 0, then all 
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points of the straight line joining this point to the point ( x , t ) are on g = 0. 
This is the property of a cone. To prove this fact, note that if (y, t) satisfies 
Eq.(16), then the point A = (ay + ( 1 - a)jc, ax + ( 1 - a )t ) where a is an 
arbitrary number, also satisfies Eq.(16). But the point A lies on the straight 
line joining ( x , t) and (y, t) in the four dimensional space. 

The wave operator in Eq.(15) has constant coefficients. Therefore, the charac- 
teristic cone and the characteristic conoid at any point coincide. Equation 
(16), thus, describes the characteristic conoid at the point ( x , t ) . Remember 
that ( x , t) is fixed in our analysis and this point specifies the position of the 
vertex of the cone and the conoid. What is the surface g = 0? Since it is 
impossible to visualize a surface in four dimensions, we must keep some 
variables fixed. Let us fix the source time variable t in addition to (jc, t) . 
Then, the surface g - 0, in variable y, is a sphere in three dimensional space 
with center at x and radius c(t-x) for all time x <t. The radius of this 
sphere becomes zero at x = t . As the source time variable x varies from - °° 
to t, the radius of the sphere shrinks at the speed of sound c from °° to zero. 
This is why the characteristic cone g = 0 of the wave equation is called the 
collapsing sphere when viewed dynamically as explained here. We will use 
this terminology in the present paper. 

From the theory of hyperbolic partial differential equations [21,22], the char- 
acteristic conoid at the point (jc, t) determines the domain of dependence of 
this point on the initial and the boundary data. We can illustrate this schemat- 
ically for the wave equation with constant coefficients in three space vari- 
ables. The equation of the characteristic conoid is given by Eq.(16). Figure 1 
shows this cone intersecting the three dimensional space, shown as planes, at 
various source time x < t. The domain of dependence of the point (jc, t ) is 
the circular region inside the cone where it intersects the three dimensional 
space at the given source time. Any data specified on this circular region 
influences the solution of the wave equation at (jc, t) . In the problems consid- 
ered here, the data are specified on a moving surface, e.g. a propeller blade, 
and our wave equation has three space variables. How do we determine the 
domain of the dependence of the space-time point (jc, t) in this case? 
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The surface / = 0 on which aerodynamic, kinematic or Kirchhoff data are 
specified, moves and possibly deforms in space. Since we have reserved the 
variables (x, t) for the observer, we now use the source variables for this sur- 
face, that is we write f(y, t) = 0. The influence surface, or the E- surface 
associated with a given observer space-time variables ( x , t) is obtained by 
the intersection the characteristic cone (conoid) of this point and the surface 
f(y, t) = 0. The mathematical description of the E- surface is obtained by 
eliminating t between / = 0 and g = 0 to get Eq.(14). The E- surface is 
easily visualizable and constructible in space. Since ( jc , t) is fixed, 
g(y, t 0 ; x, t) = 0 is a sphere for a fixed value of source time x Q . The sur- 
face f(y, Tq) = 0 describes the frozen position of the moving body or the 
Kirchhoff surface at the source time Tq . We vary source time from - °° to 

T = t and the locus of the curves of intersection of / = 0 and g = 0 forms a 
surface in space which is the E - surface for the space-time observer variables 
(x, t) . The curve of the intersection of the surfaces / = 0 and g = 0 is called 
the T -curve in our previous publications on rotating blade noise [15,16,23]. 
Figure 2 shows the collapsing sphere at a fixed source time intersecting a 
helicopter blade which is assumed to be flat. Figure 3 shows the construction 
of the E - surface for only one flat blade of a helicopter rotor in forward 
flight. The E - surface is obtained dynamically by finding the T -curve for 
several source times. Note that the E- surface is swept back for the rectangu- 
lar blade planform and lies in the rotor plane. It is flat, i.e. its mean curvature 
is zero, because we assumed a flat blade. In general, the surface /( x, t) = 0 
is not flat and the mean curvature H F of its E - surface is not zero. In the next 

section we will derive the mathematical expression for H F . 

4.0 The Derivation of the Mean Curvature of the E - Surface 

We use classical differential geometry [12-14] to derive the required mean 
curvature of the E- surface. It is crucial to select appropriate variables in 
deriving this result for two reasons. First, as will be seen below, the algebraic 
manipulations are complicated and the selection of appropriate variables 
makes these manipulations manageable. Second, we want to give geometric 
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interpretation to many of the terms in the analytic expression of the mean cur- 
vature. This requirement actually allows us to write the expression for the 
mean curvature in the most compact form. We mention here that, since the 
equation describing the £- surface is F(y\x,t ) = 0, theoretically one 
should be able to get the mean curvature H f from this equation. Such an 

approach, however, appears impractical because it does not satisfy the above 
two requirements. First, the algebraic manipulations are intractable and the 
final expression for the mean curvature is very long and complicated. This 
complexity has undesirable impact on computer code development and effi- 
ciency. Second, we get many terms whose geometric interpretation is unclear. 
We, thus, looked for variables which would make the derivation of the mean 
curvature more direct and which would satisfy the above requirements. 

4.1 The Selection of Variables 

We describe the surface f(y, t) = 0 by the time dependent position vector 
12 12 

y(u , u , x). Here, the variables (u , u ) are the Gaussian coordinates on 
the surface f(y, t) = 0 at the time t. The £- surface is described by the 
position vector: 

12 12 12 
y '(u , u ; x, t) = y[u , u , \\t(u , u ; x, 7)] (eq 17) 

1 2 

where the retarded or the emission time function \j f(u ,u ; x, t) is the solu- 
tion of the equation 


g = V -, + - 


x-y{u l , u 2 , \|/) 


= 0 


(EQ 18) 


Note that y, y ' and \| f are functions whose dependence on variables are 
clearly specified. Here the function \\f gives the emission time of the point 

1 2 

with Gaussian coordinates (u , u ) on the surface / = 0. Therefore, y ' lies 
on the £- surface. We have implicitly defined the dependence of the position 
vector y ' on y through Eqs.(17) and (18). As will be seen below, this is suf- 
ficient to give us the result that we are looking for. We remind the readers that 
in all the analysis in this section, the observer variables (x, t) are kept fixed. 
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Let R(u , u ) be a surface in the three dimensional space. Denote the coeffi- 
cients of the first and the second fundamental forms of this surface by 

g n and b , respectively. We define the matrices G and G 1 as follows: 


G = 


^11 ^12 

, G-‘ = 

11 12 
8 8 

/21 g 2Z 


21 22 
8 J 


(EQ 19) 


where G is the inverse of G. The mean curvature of the surface is given by 
the relation: 


H = \b\ = \(b\ + b 2 2 ) 


(EQ 20) 


l IK 

where bj = g . We have used the summation convention on repeated 
index here. Let g ^ = det (G) , then the element of the surface area is given 
by 


dS = Jg(2^du^du‘ 


(EQ 21) 


The above results are in books on classical differential geometry [12-14]. We 
have given a summary of the useful differential geometric results in reference 
[16]. Below, we use Eq.(20) to derive the mean curvature of the £- surface. 


Remark: We use primed symbols to denote the geometric parameters of the 
£ - surface whereas the unprimed symbols are reserved for the parameters on 
the surface / = 0 . 


4.2 Coefficients of the First Fundamental Form of the £- Surface 

Let (y^, y^) and (y ' , y ^ ) be the natural basis vectors for the surfaces 
/ = 0 and the £- surface F = 0 . These are defined by the relations: 
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From Eq.(17), we find 


v = v . + 

dx 


dy d\\ f 


du 


From Eq.(18), we get 

dy _ !i 

du‘~ 


(EQ 23) 


(EQ 24) 


where r i = r.y i , r = \x-y\/r and M r - M.r. Here M is the local 
Mach number of the surface f = 0 given by the following relation 


M = ij* 
c ox 


(EQ 25) 


Note that this vector is unambiguously defined even for a deforming surface. 
Using Eqs.(24) and (25) in Eq.(23), we get 


y = + 


r i 


1 -M 


M 


(EQ 26) 


From this relation, we construct the coefficients of the first fundamental form 
of the E - surface in terms of known quantities on the surface / = 0 as fol- 
lows: 


M.r. + M.r. M 2 f.r. 

, „ i j j i , i j 

6 ij 6 ij 1 _ M 2 

r (1 -M r ) 


(EQ 27) 


where g.. = y y is the coefficients of the first fundamental form of the 
surface / = 0, M- - M.y- is the covarient component of M with respect to 
the natural basis vectors on the surface / = 0, and M = \M\ . Note that g' ■■ 

u 
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depends on the geometric parameters of the surface / = 0, the kinematic 
parameters as well as the radiation direction f . . 


From Eq.(27), we find the determinant of the coefficients of the first funda- 
mental form of the £ - surface . The derivation of this result involves much 
algebraic manipulations. The final result, however, is simple and it is 



A A 2 
1 ~M r ) 8 (2) 


(EQ 28) 


where A is given by Eq.(12) and g is the determinant of the first funda- 

(2) 

mental form of the surface / = 0 . Now we use this equation to relate the 
element of surface areas d.L and dS of the £ - surface and the surface 
/ = 0 , respectively. This relation is obtained from the determinants of the 
coefficients of the first fundamental forms on these surfaces and using 
Eq.(21). It is 




(2) 


du^du^ = 


A 

II -Af 


■f* 


(2) 


du^ du^ 


A 

1 -M r 


dS 


(EQ 29) 


This very useful result was originally obtained by the first author using a geo- 
metrical method [23]. 


We can now find the elements g' lJ of the inverse of the matrix 
G ' = [g' d .These are defined as follows 

U 


g 


11 


8 22 
§ '( 2 ) 


.12 „’ 21 

g = g 




S'll 

(EQ 30) 

8 (2) 


We have given g’-- and g’^ above in Eqs.(27) and (28). We introduce two 
new vectors tangent to the surface / = 0 : 

y = n x r , X = n x M (eq 3i) 
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in order to write g' IJ in contravariant tensor form. We note that the contravar- 
iant components of the vectors y and X are given by the relations 



respectively. We can then write g' lJ as follows: 




(1 ~M r ) ij 1 -M 


A 


„ r , i~j M i j 

-g + — 7- (Y A + y J X ) + — y Y 
A A 


(EQ 32) 


(EQ 33) 


(EQ 34) 


We will use this result in the next subsection. 


Remark: We clarify an important point about notation here. What is the 
meaning of the parameters of the surface / = 0 on the £ - surface in the 
above equations? The L- surface is given by F(y,x,t ) = 0 where the 
observer variables (x, t) are fixed. An emission time is associated with each 
point y on this surface. This emission time is given by t -\x - y\/c . At this 
instant, the surface / = 0 intersects the £ - surface along a curve which we 

called the F -curve (see Figs. 2 and 3). All parameters of the surface / = 0 
at the point y on the £ - surface are calculated at the emission time of the 
point y. This means that a parameter such as M f , in Eq.(29) is evaluated as 
follows: 


M r {y\ x, t)=M^y,t- \— ^ y ^ j.r (EQ35) 

where M and r = (x - y)/\x - y\ are calculated from the surface / = 0 

but precisely on the Y -curve it shares with the £- surface . Note that all the 
variables (y, x, t ) on the right of this equation are independent variables and 
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there is no ambiguity in the interpretation of the dependent variables as we 
have written them here. As explained below, differentiation of the dependent 

1 2 

variables with respect to variables ( u , u ) , however, needs more care. 


4.3 Coefficients of the Second Fundamental Form of the E- Surface 


The coefficients of the second fundamental forms of the surface / = 0 and 
the E- surface are defined by the relations: 


b 

ij 


n 


d y 

dudiJ 


n.y 

ij 


b ' 
ij 



dudu J 



(EQ 36) 


where N is the unit normal to the E - surface which is given by 


N = 


VF 

|VF| 


n-M r 

n 

A 


(EQ 37) 


The subscripts of y and y ' denote differentiation with respect to variables 
1 2 

(u , u ) . As will be shown below, the derivation of the analytic expression 
for b requires considerable amount of algebraic manipulations. The reason 

1 2 

is the dependence of the emission time on the variables (u , u ) . Thus, for 
any independent variable defined on the E- surface , e.g. n above, we have 


dn 3,1212 3,12 dn 3 \\f 

— : = — : n(u , u , \| f{u , u ; x, t)) = — : n(u , u , \|/) + ^ A (eq 38) 

'x l -\ l -\ l C/T 'x l 

ou ^ ou ou ou 


where the emission time function \| f is given by Eq.(18). Here, on the right of 

the second equality sign, we have taken the unit normal n as a function of 

12 12 12 12 
variables ( u ,u , x). Note that n(u ,u,t) and n(u , u , \| f(u ,u\x, t )) 

are two different functions- the former is defined on the surface / = 0 and 
the latter on the E - surface . Strictly speaking, we should use different sym- 
bols for these functions but this is not done in practice because it leads to pro- 
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fusion of symbols which can be annoying. We should warn the readers that 
this practice, common as it is, can lead to confusion at times in wave propa- 
gation problems. For this reason, we have made an exception to the usual 
practice for one symbol that appears very often in our work. We introduced 
the symbol x\f for the retarded time function in place of which we could have 
used the same symbol as the source time x. We have reserved t for an inde- 
pendent variable. 


To reduce the risk of confusion in interpreting a partial differentiation with 
1 2 

respect to (u , u ) , we specifically indicate this on the £- surface as follows 


d/ du 


£' 


This means that the independent variable is a function of 


( u , u ; x, t ) . However, we use the notation d/du for partial differentiation 

1 2 

of an independent variable which is a function of (u , u , x) on the surface 
/ = 0. There is no confusion for the partial derivative d/dx which applies 
only to the independent variables on / = 0 . Note that on the £ - surface , 
differentiation with respect to the source time is meaningless since the inde- 
pendent variables on this surface are not functions of t . 


By taking the derivative with respect to the variable u ] of both sides of Eq. 
(23) and after much algebraic manipulations and simplifications, we get 



s 'a M 

’v r( 1 - M r ) 


+ 


1 

1 -M 


f 


dM „ dM 
r • — : + r •- 


\ 


J- 


\ 'd 


u 


du 


J 


J 


+ 


M 

1 -M r 



+ b-. cos0) 
l J 


+ 


f.f.M 
l J 


+ 


M 


c(l-M r ) 2 (1 -M r ) 2 


1 , dM , dM'' 
rr : + r r : 

J -S l 


V 


du J 


du 


Mr.f.M f.f.M 

^ J _i_ ^ J 

c(l-M r ) 3 r( 1 - M r ) 3 


(EQ 39) 
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k 

Here T ■ ■ is the Christoffel symbol of the second kind, 0 is the angle between 

n and the radiation direction r, and M r = r M . We have summed some 
terms on the right sides to get g' ■■ in the numerator of the second term 
because in obtaining the mean curvature of the £ - surface , we must multiply 

b' •• with g' l] as will be seen later. Therefore, this summation reduces some 
algebraic manipulations. 


4.4 The Mean Curvature of the 2-Surface for a Rigid Moving Surface 


So far we have not assumed that the surface / = 0 is rigid. All the above 
results are valid for a deformable surface. Now, we will make the assumption 
of rigidity of the surface. In applications that we have in mind, the surface 
/ = 0 is always rigid, e.g., the surface of a propeller or rotor blade. The 
Mach vector distribution of the surface can be written as 
M = (v 0 (T) + coxr|)/c where v 0 (t) is the velocity of a reference point on 

the surface, r\ is the position vector of a point on the surface from the refer- 
ence point, and to is the angular velocity of the surface. We note that when t 
is fixed, we can assume that y and r\ coincide and, thus, we get 

dM/du = to x y/c. Using the definition of the second fundamental form, 
Eq. (36), we obtain from Eqs. (37) and (39): 



b u , f i q j + f j q i 
A A(1 -M r ) 


+ 


r.r M 
i j n 


cA(l-M r Y 


M g' ■■ 
» 6 u 

rA 


+ 


r.r .M 
i J n 


rA( 1 - M A 


where we have defined 


(EQ 40) 


1 

q = -ii x co 


(EQ 41) 
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and A is given by Eq. (12). We are now ready to find the mean curvature of 
the £ - surface for a rigid moving surface. 

From the definition of the mean curvature, Eq. (20), using primed variables 
for the £ - surface , we obtain the following expression for H F : 


_ M n( sin 2 0A (l 

H f TaI 1 _ TTT J + .3 H f 


2A 


A' 


^ ^ r( ^ ^"1-1 y2~2 N 

+ — r-g + KjX y + k 2 A y 

A V 


+ 


sin 2 0(T ■ 2 

T \- M n + M K 

2A vc 


+ A(Y q)(k r) 

A 


Here we have used the following definitions of the symbols: 


A 

o 1/2 

= (1 +M -2M cos6) 

v n n ' 

COS0 

-nr 

r 

- (x - y)/r unit radiation vector 

n 

= V/ (unit outward normal to/= 0) 

M 

= dy/cdx , M = \M\ , M n = M n 

H t 

= local mean curvature of/= 0 

Ky and 

k 2 = local principal curvatures of/= 0 

K y 

= local normal curvature off= 0 along y 

y 

= nxr 

A 

= n xM 


(EQ 42) 


q = n x to / c 

CO = angular velocity of/= 0 

_1 ^2 ~1 ~2 

(y , y ) and (A. , X ) are the components of y and A (see Eq.(31)) in princi- 
pal directions with respect to unit basis vectors 

Equation (42) is the main result of this paper. This result was derived and 
announced by the first author in 1996 [16] but the derivation has not been 
published. The second author independently verified the final result through 
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a slightly different approach. The extension to the deformable case was also 
carried out by the second author which will now be presented. 

4.5 The Mean Curvature of the Z-Surface for a Deformable Moving 
Surface 

If we assume the surface / = 0 is deformable, then we get a slightly more 
complicated expression for the mean curvature of the £ - surface . In the deri- 
vation of the mean curvature formula, Eq.(42), the rigidity hypothesis was 
first introduced during the computation of the second fundamental form, 
Eq.(40). In order to avoid this restriction, we begin by rewriting Eq.(39) as 
follows: 


, d\\fdM 3\t/ dM 

rtf -w } ~\,.J 


+ 


M 


r( 1 - M ) 


du du J du J du 
2 chj/ chj / \ 

dudu^ 


-gy + c 


Here T is a linear operator defined by 


r(V) = v+ r^r M 


(EQ 43) 


(EQ 44) 


for a given vector V. Next we compute b' -j = y' -j ■ N . To do this, note that 
we can show T(V) • N = V n/A . From Eq. (43), we get 


b 'a = x 


1 


b[ i + 1 - M. 


T c)M „ dM 
r - — : • n + r ■ — - • n 


l ~ j i 


du J 


du 


+ 


1 i* j 


dM M J 

n h 


c(l-M r ) 


2 dx 


V 


Zij + 


T i r j 


M 


(1 


(EQ 45) 


Considerable algebra is now involved to get the mean curvature of the 
£- surface . Taking advantage of the symmetry of the hrst and the second 
fundamental forms of the surface / = 0 and the orthogonality of the vectors 
r and y yield: 
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(EQ 46) 



Here the operator L is based on the coefficients of the second fundamental 
form defined as L(U, V ) = b--UW J and the vector r is the projection of 
the unit radiation vector r on the local tangent plane to the surface / = 0 . 

/V 

Note that r f - sin 0 . 

For a deformable surface, a new problem appears which must be addressed 
carefully. Note that for the case of a rigid surface, all quantities on the right 
side of Eq. (42), or equivalently Eq. (46), all surface coordinates are indepen- 
dent and all the velocities are unambiguously defined. In the case of a 
deformable surface, we must show that if two surface Gaussian coordinate 
systems slide with respect to each other, then the mean curvature of the 
E - surface is independent of the tangential Mach numbers of the two coordi- 
nate systems. The terms of the mean curvature expression for the E - surface 
of a deformable surface can be written in such a way that the tangential Mach 
number of the coordinate system is taken out. In order to accomplish this, we 
decompose M into its normal and tangential components on the surface 
/ = 0 denoted M n and M f , respectively. We then separate the mean curva- 
ture expression into pieces that involve these normal and tangential compo- 
nents separately. We also remove the explicit dependence on the vectors y 
and X by using the identity 

L(M f r t ) + L(y,X) = 2 H f M f -r t (EQ47) 

This allows the mean curvature of the E - surface to be written as 
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(EQ 48) 


M Hr 

H p t - + —r~ + 

* rA A 


2 A' 


2 f M 2 

sin 0^— ^ - M n K r + M n 


1 , . dn ,,2 
+ - m ,-^ + m , k m 


-2(l-M„cos e/^ + L(f,,M ( ) 


where k and are the normal curvatures of the surface / = 0 in the 
directions of r f and M f , respectively. We have used the fact that 
2 

L(V,V) = \V\ Ky, where Ky is the normal curvature of / = 0 in the 
direction of the tangential vector V. 


In Eq. (48), everything is clearly independent of frame of reference except, 

1 ~ dw 

apparently, for the expressions -r ■ t— + L(r , M.) and 

c 1 OT IT 

M n + -M f ■ ^ k m .It turns out that these expressions are also indepen- 
dent of frame of reference. In fact, it can be shown that 

-Jt ■ |! + L ^r ■ V / T (EQ 49) 

and 



1 ,, dn ,,2 
+ - M . • + M . k 




t M 


\ _ f 1 

2-^xx + c \ 


(EQ 50) 


Although Eqs. (49) and (50) show that the mean curvature expression Eq (48) 
is independent of the coordinate system chosen, Eq. (48) as it stands might be 
more suitable for use in practice. 

Now that we have derived the analytic expression for the mean curvature of 
the £- surface, we can write the solution of the FW-H equation and the 
Kirchhoff formula for a supersonically moving surface in a very compact 
way using Eqs. (11) and (13) in the solution of the wave Eqs. (4) and (8). We 


The Mean Curvature of the Influence Surface of Wave Equation With Sources on a Moving Surface December 8, 1999 


24 



feel that by relating the mean curvature of the E - surface to the geometric 
and kinematic parameters of the surface / = 0 , we have derived useful for- 
mulas for applications in aeroacoustics. 


5.0 Discussion and Concluding Remarks 

The mean curvature of the E - surface is a very complicated expression and 
there is little hope that for moving surfaces of engineering interest one can 
obtain an analytic result given the geometry and the history of the motion of 
the moving surface / = 0 . A simple but not sufficient test of validity of H F 
is that for a stationary surface, one should get H F = Hj. This can be easily 
established from Eq. (42). However, numerical test of this result is possible 
for a given smooth surface in rectilinear and rotating motion. Note that the 
mean curvature of the collapsing sphere is 1/r. As seen from Eq. (42), only 
the first term of the mean curvature of the E - surface depends on the mean 
curvature of the collapsing sphere. 

One of the main reasons for studying the mean curvature of the E - surface 
has been to explain the complexity and the appearance of normal curvatures 
of the surface f - 0 in many directions in Formulation 3 of Farassat used 
for supersonic propeller noise prediction [11]. It is now clear that Formula- 
tion 3 can be written in a much more compact form using the mean curvature 
of the E - surface . The procedure for the derivation of this compact formula 
has been outlined in Section 2 of this paper. Such a formula will simplify 
coding for noise prediction by calculating H F in a separate subroutine. We 
have made further progress in finding a very simple solution of Eq. (10) 
which does not depend on the mean curvature of the E- surface [24]. We 
feel, however, that the solution involving the mean curvature discussed here 
is fundamental and it will lead to further understanding of noise generation 
from moving surfaces. 

There is one important point that must be mentioned in connection with the 
solution of the wave equation of Eq. (10). When the collapsing sphere leaves 
the surface f = 0 at a point where M n = 1 , then one gets A = 0 in the 

denominator of the integrand of the solution of Eq. (10) given by Eq. (13). 
This can lead to infinite singularity in the acoustic pressure p' . Based on a 
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recent study by Farassat, Brentner and Dunn [24], it can be shown that in Eq. 
(13), there is a cancellation of the highest order of 1/A singularity between 
the first and the second terms in the integrands on the right side of the equa- 
tion. Because of the particular structure of q(y, t) in the governing equation 

of the Kirchhoff formula [24], the acoustic pressure calculated from the 
Kirchhoff formula is finite for all time. Similar situation exists for the FW-H 
equation if we add the surface terms coming from the quadrupole source term 
to the thickness and loading terms [24]. However, noise calculation based on 
the thickness and loading terms alone for supersonic blades with blunt lead- 
ing edge can result in infinite acoustic pressure. 
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Figure 1. Schematic diagram of the characteristic cone (conoid) of the 
observer space time variable (x, t ) and its domain of dependence (DOD). 
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Figure 2. The collapsing sphere intersecting a rotating blade at the source 
time T . Note that the collapsing sphere is the intersection of the characteristic 
cone (conoid) of ( x , t ) with the three dimensional space at the time T . 



Figure 3. The construction of the £- surface for a flat rectangular rotor 
blade. The rotor is moving forward at the velocity V F in such a way that all 

points on the rotor travel at subsonic speed relative to the speed of sound in 
undisturbed medium. The T - surface will be more complicated and even in 
several pieces if part or all the blade surface travel at supersonic speed 
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